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Abstract. We evaluate the complex frequencies of the normal modes for the 
massive charged scalar field perturbations around a Reissner-Nordstrom black hole 
surrounded by a static and spherically symmetric quintessence using third order WKB 
approximation approach. Due to the presence of quintessence, quasinormal frequencies 
damp more slowly. We studied the variation of quasinormal frequencies with charge 
of the black bole, mass and charge of perturbating scalar field and the quintessential 
state parameter. 
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1. Introduction 

It is now well understood that different kinds of perturbations in the geometry of a 
black hole can excite certain combination of its characteristic complex frequencies of 
the normal mode, called quasinormal modes(QNMs), whose real part represents the 
ring down frequency and imaginary part, the decay time and they are independent 
of the initial perturbation, but depends only on the parameters of black hole. The 
concept of QNM is put forward by Vishweshwara[2] and since then the quasinormal 
spectrum of black holes has been extensively studied for a great verity of black hole 
back grounds and perturbation fields because of its astrophysical and other theoretical 
interests. Astrophysical interests are associated with their relevance in gravitational 
wave analysis[3]. The QNMs of black hole are expected to be detected by the future 
gravitational wave detectors such as LISA[15], and give an opportunity to explore the 
properties of black holes directly. Apart from the observational interest of detection of 
quasinormal ringing by gravitational wave detectors, study of QNM found significance 
in AdS/CFT correspondence [17], black hole area quantization and Loop Quantum 
Gravity [16J. 

The idea that our universe is in a phase of accelerated expansion rather than 
holding steady, is strongly supported by a set of recent interlocking observations such as 
studies of distant super novae pQ, type la supernova[3] and cosmic microwave background 
radiation anisotropyp)], indicating the presence of some mysterious form of repulsive 
gravity called dark energy. In order to explain the nature of dark energy several models 
were proposed(for a recent review see [6]). The simplest option for this dark energy is 
Einstein's cosmological constant [7] with a constant equation of states(EOSs) e = — 1 but 
it needs extreme fine tuning to account for the observations. The second is the dynamical 
scalar field models like quintessence [8j, k-essence[9] and phantom[10], in which the EOSs 
of dark energy changes with time. Among them, the quintessence is the most natural 
model, and can give the EOSs with — 1 < e < —1/3. 

Because of relevance of the two topics, some studies of QNMs of black hole were 
started in the presence of quintessence, after Kiselev[18] derived the exact solution 
of Einstein equation for quintessential matter surrounding a black hole. QNMs of 
Schwarzschild black hole is studied for scalar [TT], gravitational[12], electromagnetic [13] 
and massless Dirac fieldjH] perturbations in the presence of quintessence. Our work 
extend the studies to charged black hole space times. Charged scalar QNMs of 
Reissner-Nordstrom black hole were studied in[23[ [23]. In this paper we considered the 
perturbation of Reissner-Nordstrom black hole encircled by quintessence for massive 
charged scalar field perturbations and study the dependence of quasinormal spectrum 
upon different parameters of the problem such as charge of the black bole, mass and 
charge of perturbating scalar field and the quintessential state parameter. 
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2. Massive charged scalar field around a Reissner-Nordstrom black hole 
surrounded by quintessence 

Kiselev derived a static spherically symmetric exact solution of Einstein equations for 
quintessential matter surrounding a charged black hole under the condition of additivity 
and linearity in energy momentum tensor [18J. The metric can be expressed in the form 



ds 2 = -f(r)dt 2 + f(r)- l dr 2 + rW 



(1) 
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with M and Q mass and charge of the black hole, e is the quintessential state 



parameter and c is the normalization factor which depends on p q 
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of quintessence. The Klein-Gordon equation describing the evolution of massive charged 
scalar perturbation field outside a charged black hole is given by [19] 

$ ;/ ^ - ieA^(2^. v - ieA v <$>) - ieA^g^Q = m 2 $, (2) 

where is the electromagnetic potential and e is the charge of the scalar field. 
Representing the charged scalar field in to spherical harmonics 
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the wave equation becomes 
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where we used the coordinate transformation defined by 
dr 



n - e 2 A 2 r] = 0, 
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The electromagnetic potential A t = C — ® , where C is a constant. We define 
T] = e~ ieCt y, 



(4) 



(5) 
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to avoid the physically unimportant quantity C. where $ is the radial part of 
perturbation variable taken to have time dependence e~ tu}t . Now the radial perturbation 
equation can be written as 



+ 6^ = 
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where G 
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-V 2 and V is the scattering potential, which is a function of frequency 



of perturbation, u and is given by 



V = f(r) 
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3. Evaluation of Quasinormal modes 

We use the third order WKB approximation method to determine the complex normal 
mode frequencies of black hole, a semi analytic method originally developed by Schutz 
and Will [20] for the lowest order. Later Iyer and Will [21 J carried this approach to third 
WKB order and Konoplyap2] to sixth order to get more accurate results. 

It gives a simple condition which will lead to discreet, complex values for the normal 
mode frequencies. 

A %-A(n)-n(n)=n + i, (9) 
26o 2 

where A and Q are higher order terms given by 

A <"> = <*M (f ) (i + Q ' 2 ) - m (S) 2 ( 7 + 6te2 )i- < 10 > 



n 




(11) 



where 



1 ... r 0,1,2, Re(E)yO 

a = n-\ — with n = \ . . 

2 1 -1,-2,-3, Im(E) -< 

Qq^ denotes the n th derivatives of G evaluated at r , the value of r at which V 
attains maximum. Here a complexity arises from the fact that the potential V is a 
function of frequency u. This makes difficult to calculate the value of Tq- We make 
use of the procedure suggested by Konoplya[23] to find r by fixing all the parameters 
other than u for which the V depends and then find the value of r at which V attains 
maximum as a numerical function of u. Substitute for ro and then find the value of u 
which satisfies the condition (Q by trial and error way. 



4. Results and Discussion 



We take all the values of parameters in black hole mass(M) units. It is a general 
experience that the imaginary part of quasinormal frequencies decreases with increase 
in mode number, n, means that the quasinormal frequencies with lower mode number 
will decay slowly and are relevant to the description fields around black hole. So we 
consider frequencies for lower mode number for our study. The dependence of real and 
imaginary parts of u on the charge of the black hole, Q is plotted in Figur^H for fixed 
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I = 3, n = 0, e = 0.1, m = 0.1 and for different values of e with c = 0.001. The case with 
absence of quintessence is also plotted (dotted line). 
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Figure 1. Reuj and Imui as a function of Q for I = 3, n — 0, e = 0.1, m = 0.1 and 
for different values of e with c = 0.001. The dotted line represents the no quintessence 
case(c=0). 

The quasinormal frequencies for scalar field in charged black hole is influenced by 
quintessence. The magnitudes of real and imaginary parts of uj is lower in the presence of 
quintessential field. This implies due to the presence of quintessence, the quasinormal 
frequencies for scalar field in RN black hole damps more slowly. As in the absence 
of quintessence [23l El], Re(u) increases monotonically with the increase in Q while the 
magnitude of Im(u) first decreases, falling to a minimum around Q = 0.8 and thereafter 
increases sharply. For larger values of Q, the modification with quintessence decreases. 

FigureS] shows the variation of Reu and Imu with quintessential state parameter 
e with c = 0.001 for fixed I = 3, n = 0, e = 0, m = 0.1. As the value of e increases Reuo 
increases slowly but the magnitude of Imu increases more fast, means damping is less 
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for lower values of e. 
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Figure 2. Variation of Reoj and Imuj with e for I = 3, n = 0, e = 0, m = 0.1 

In Figur^S] Reu and Imuj are plotted as functions of e with I = 3,n = 0, m = 0.1 
for Q = 0.1,0.2,0.3 and different values of e. Dotted line represents absence of 
quintessence(c=0). The variation is almost linear. In the presence of quintessence, 
the magnitude of Reu and Imu increases with e imitating the no quintessence case but 
with lower values of Reuj and \Imuj\. The modification with e becomes larger for higher 
values of Q . 



s 





Figure 3. Variation of Retu and Imu with e for I = 3, n = 0, m = 0.1, Q = 0.1, 0.2, 0.3 
and different values of e. Dotted line is for c = 0. 



Finally, we study the role of mass of scalar field on quasinormal frequencies. 
Quasinormal modes occur only when the peak value of the potential V(r = r max ) is 
larger than m 2 and uj 2 of the field is smaller than this peak value[25j. This means 
that there exists a maximum value for mass, m max beyond which there will be no 
quasinormal modes. m max can be estimated from the condition for the existence of 
quasinormal modes, 



(12) 
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The values of m max obtained for different values of e is tabulated in Tabl^U In the 
presence of quintessence m max decreases because quintessence lowers the height of the 
potential barrier as shown in Figure^] and when e = —1, it has the lowest value. 
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0.001 


-1/3 


0.88516 
0.88366 


0.001 
0.001 


-2/3 
-1 


0.87815 
0.85864 



Table 1. The limit of mass of scalar field, m max for the existence of quasinormal 
frequencies with I = 3, e = 0, Q = 0.1 

The values of Reuo and Imuo evaluated using WKB method with the variation of 
mass is plotted in Figured 




Figure 4. Variation of Reuj and Imtu with e for I = 3, n = 0, m = 0.1, Q = 0.1, 0.2, 0.3 
and different values of e. Dotted line is for c = 0. 




Figure 5. Effective potential I = 3,n = 0, m = 0.1, Q = 0.1,0.2,0.3 and different 
values of e. Dotted line is for c = 0. 
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Reu increases with increase in mass while |imco>| decreases with mass, which 
indicates that QNMs of fields with lower mass damp slowly. But they behave abnormally 
near m max . This is due to the fact that, for larger mass of the field the potential looses its 
barrier shape by broadening the potential peak as shown in Figur^5] and WKB method 
gives inaccurate results. Lower modes show less abnormality showing that WKB method 
is more accurate for fundamental modes. 

Another thing we noticed is that this abnormal behavior is lower in the presence of 
quintessence and even when e = —1 we can get satisfactory curve because the peak of 
effective potential broaden much less in the presence of quintessence comparing with the 
normal case(c=0) as understood from FigureSland quintessence helps to retain barrier 
shape and give more accurate results at larger mass range. 
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